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ABSTRACT 

We present Keck/HIRES spectra of 3 globular clusters in the outer halo of M31, 
at projected distances beyond « 80 kpc from M31. The measured recession velocities 
for all 3 globular clusters confirm their association with the globular cluster system of 
M31. We find evidence for a declining velocity dispersion with radius for the globular 
cluster system. Their measured internal velocity dispersions, derived virial masses 
and mass-to-light ratios are consistent with those for the bulk of the M31 globular 
cluster system. We derive old ages and metallicities which indicate that all 3 belong 
to the metal-poor halo globular cluster subpopulation. We find indications that the 
radial gradient of the mean metallicity of the globular cluster system interior to 50 kpc 
flattens in the outer regions, however it is still more metal-poor than the corresponding 
field stars at the same (projected) radius. 
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1 INTRODUCTION 

One of the great unanswered questions in astrophysics is the 
formation and evolution of galaxies in the Universe. Since 
M31 is our closest spiral galaxy, located at a distance of 
~ 780 kpc (Holland 1998), it is considered an ideal astro- 
physical laboratory to test the different ideas about galaxy 
formation and evolution. 

In particular, the outer regions of galaxies are expected 
to hold a wealth of information about the way in which 
galaxies are assembled. For example, in order to explain how 
the Milky Way was formed, Eggen, Lynden-Bell & Sandage 
(1962) and Searle & Zinn (1978) studied the kinematics and 
metallicity of stars and globular clusters in the bulge/halo 
of the Milky Way. While the former suggested the Galactic 
halo formed during a rapid (~ 10 8 years) monolithic col- 
lapse of a gaseous protocloud, the latter suggested formation 
via the accretion and merging of independent protogalactic 
fragments over a longer dynamical timescale (10 9 years). Im- 
portantly, Searle & Zinn's (1978) seminal paper on Galactic 
globular clusters constitutes the basis for modern ideas re- 
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lating to the formation of cosmological large scale structures 
in the universe (e.g. White & Rees 1978). 

Perhaps the main difference between the globular clus- 
ter systems of M31 and the Milky Way is that the former 
is more populous with an estimated total of ~450 members 
(Barmby et al. 2000; Perrett et al. 2002; Huxor et al. 2008; 
Caldwell et al. 2009) and may contain a significant popula- 
tion of intermediate age (3-6 Gyr) globular clusters. Other- 
wise they share many similarities: a Gaussian-like luminosity 
function that peaks around Mv = -7.62 mag with a disper- 
sion of 1.06 mag (Barmby, Huchra & Brodie 2001); two sub- 
populations of mean metallicity [Fe/H] = -1.57 and -0.61 
associated with the halo and bulge, respectively (Barmby, 
Huchra & Brodie 2001; Forbes, Brodie & Larsen 2001); and 
similar structural properties (Barmby et al. 2007). 

A relevant recent discovery is that of a metal-poor stel- 
lar halo in M31 (Kalirai et al. 2006; Chapman et al. 2006). At 
projected radii beyond 60 kpc Kalirai et al. (2006) estimate 
a mean halo metallicity of [Fe/H] = -1.26, which decreases 
to -1.48 if an [a/Fe] abundance ratio similar to the Milky 
Way's stellar halo is assumed. Beyond a projected radius of 
70 kpc, outer halo globular clusters in M31 have also been 
discovered recently but remain rare, these include 1 globu- 
lar cluster reported in Martin et al. (2006), 2 in Mackey et 
al. (2007), 2 in Galleti et al. (2007) and 1 in Huxor et al. 
(2008). 

In this Letter, we present Keck spectra for some of the 
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Table 1. Properties for the 3 outer globular clusters. 



Properties* 


GCM06 


GC5 


GC10 


(1) 


(2) 


(3) 


(4) 


R.A. [J2000] 


00:50:42.5 


00:35:59.7 


01:07:26.4 


Dec. [J2000] 


32:54:59.6 


35:41:03.6 


35:46:49.7 


_E B -v[mag] 


0.08 


0.08 


0.09 


rh [pc] 


2.3 ± 0.20 


6.3 ± 0.15 


4.3 ± 0.15 


M v ,o [mag] 


-8.5 ± 0.3 


-8.8 


-8.3 


Rp [kpc] 


100 


78.5 


99.9 



Notes. — (*) As given in Martin et al. (2006) and Mackey 
ct al. (2007). 



furthest projected distant (R p ) globular clusters known at 
this time in the halo of M31 — GC5 and GC10 (Mackey et al. 
2007) and GCM06 (Martin et al. 2006). Two of them (GC5 
and GC10) do not have spectroscopic information published 
to date. Thus, we present, for the first time, their kinematics, 
structural parameters, ages and spectroscopic metallicities. 
We briefly discuss our results within the context of the halo 
assembly of M31. 



2 OBSERVATIONS AND DATA REDUCTION 

The three globular clusters were selected from the imag- 
ing survey reported by Martin et al. (2006, GCM06) and 
Mackey et al. (2007, GC5 and GC10). The spectra were ob- 
served with the Keck/HIRES instrument (Vogt et al. 1994) 
on 2008 August 19. The instrument setup employed a 0'.'86- 
wide, /.'O-long slit, providing a spectral resolving power 
of R w 50000. The cross-disperser angle was set to cover 
wavelengths AA = 4020-8520 A. For each globular cluster, 
3x900-s exposures were taken in succession, without any 
intervening re-acquisition sequences, with the slit oriented 
at the parallactic angle. Exposures were offset from each 
other in the direction perpendicular to the parallactic angle 
(i.e. along the horizon). The first exposure was centred on 
the globular cluster with the others offset either side of cen- 
tre. The offsets were l'.'O for GC5 and GCM06 and 0'.'8 for 
GC10. The seeing was stable throughout the observations 
at ~078. 

Science and calibration data were reduced using the 
MAKEE (MAuna Kea Echelle Extraction) Keck Observa- 
tory HIRES Data Reduction Program, written by Tom Bar- 
low. Bias frames, flat-fields, order trace and ThAr arc images 
were taken as part of the Keck base calibrations. The final re- 
duction process comprised bias-subtraction and flat-fielding, 
using base calibrations which were previously co-added. The 
wavelength calibration was carried out with solutions ob- 
tained from the ThAr arc exposures, which provided typi- 
cal mean residuals of 0.05 A. Spectra were then cosmic ray 
cleaned, sky-subtracted and extracted into a series of helio- 
centric corrected 2D spectra with no correction to a vacuum 
wavelength scale applied. The spectra were then combined 
with inverse variance weighting to yield the final spectrum 
of each cluster. Basic properties for the 3 globular clusters 
are given in Table Q] Figure [T] displays a portion of the final 
spectra indicating the data quality. 




5160 5165 5170 5175 5180 5185 

Wavelength [A] 



Figure 1. Integrated-light Keck/HIRES spectra of the 3 globular 
clusters around the Mgfe Triplet wavelength region. The promi- 
nent absorption Mgl lines are labeled just below the spectra. The 
typical signal-to-noise per pixel is 20 at this wavelength region. 



3 DATA ANALYSIS 

The galactocentric coordinates (X,Y) and projected dis- 
tances for our 3 outer halo globular clusters were calcu- 
lated relative to an adopted M31 central position of cmooo = 
00 h 42 m 44 s .31, <5j2ooo = +41°16'09 / .'4. The position angle for 
the X-coordinate is 38° (Kent 1989), and we adopt 780 kpc 
as the distance to M31 (Holland 1998). At this distance, 1' 
corresponds to 228 pc. The projected distances calculated in 
this way are similar to those quoted in Tabled] from colour 
magnitude diagram (CMD) studies. These calculated radii 
are used throughout this paper to facilitate comparison with 
other literature work. 

The velocity dispersions (er) and heliocentric radial ve- 
locities (vq) were obtained by fitting the observed spectra 
with the convolution of suitable template spectra of giant 
stars taken with the same telescope and instrument and with 
a Gaussian velocity profile using the program pixfitgau 
(see van der Marel & Franx 1993 for more details). The 
best fitting parameters are obtained by ^-minimisation pro- 
cedures. We used the spectral region 5150-5240 A, where 
strong features permit a reliable determination of the veloc- 
ity dispersion. The velocity dispersions were also measured 
using the spectral region 5370-5420 A, which contains many 
iron lines. The results were consistent displaying differences 
of less than 1 km s . The adopted values are based on the 
5150-5240 A spectral region. The uncertainties quoted are 
only due to the random noise, that is, by taking into account 
the CCD read-out noise per pixel (~ 0.6 e~/pix) and gain 
(1.9 e"). 

The globular cluster masses were computed by using 
the Virial Theorem, being M = Wxa 2 ^/ G, where a is the 
velocity dispersion, rh is the half-light radius, G is the grav- 
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itational constant and a virial coefficient of 10 was used (see 
Forbes et al. 2008). For GC5 and GC10 the r h values were 
measured using images obtained with the Advanced Cam- 
era for Surveys on Hubble Space Telescope (Mackey et al. 
2007), while for GCM06 it was obtained using ground-based 
observations (Martin et al. 2006). 

Spectroscopic metallicities were obtained by using 
metallicity calibrations from Mgb and CH (Brodie & Huchra 
1990; Perrett et al. 2002) and for Mg2 (Brodie & Huchra 
1990; Buzzoni, Gariboldi & Mantegazza 2002) indices. In 
addition, we also employed the robust approach of Proctor, 
Forbes & Beasley (2004), which performs the simultaneous 
^-minimisation of a large number of spectral indices. This 
method leverages the entire suite of available indices to esti- 
mate stellar population parameters and, in particular, makes 
apparent those indices which deviate from the general trend 
of the data, effectively eliminating errors in calibration or the 
spectra themselves (e.g. sky-line contamination, bad pixels 
etc). 

The Single Stellar Population (SSP) models used here 
are from Thomas, Maraston & Korn (2004), which have been 
used to recover reliable ages and metallicities for both Galac- 
tic and extragalactic globular clusters (e.g. Mendel, Proctor 
& Forbes 2007; Beasley et al. 2008). However, we note that 
such an analysis is usually applied to much lower resolu- 
tion data than obtained here. Hence, absorption line indices 
were measured after broadening the spectra to the wave- 
length dependent Lick/IDS system resolution as described 
in Worthey & Ottaviani (1997). The Lick index definitions 
were taken from Worthey & Ottaviani (1997) and Trager et 
al. (1998). 
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Figure 2. Heliocentric radial velocities for M31 globular clusters 
corrected to the M31 velocity frame against the projected radius. 
Symbols are as follows: the crosses for globular clusters at R p ^ 
40 kpc, while the solid symbols are for those at large projected dis- 
tances (R p 40 kpc). The present sample is represented by solid 
triangles, while the outer globular clusters of Galleti et al. (2007) 
and Mackey et al. (2007) are represented by the solid circles. 



4 RESULTS AND DISCUSSION 

We derive a heliocentric radial velocity of —358.3 ± 1.9 km 
s" 1 for GC5, v e = -291.2 ± 2.1 km s" 1 for GC10 and 
a value of -354.7 ± 2.2 km s" 1 for GCM06. All 3 globu- 
lar clusters are consistent with being part of the M31 (vq 
= -300 km s _1 ) globular cluster system. Recently, using a 
low-resolution spectrum (R ~ 1300), Galleti et al. (2007) 
determined v Q = -312 ± 17 km s^" 1 for GCM06. This mean 
difference between the heliocentric radial velocities suggests 
the errors quoted in the lower resolution work may have been 
underestimated. 

As shown in Fig. [2] our radial velocity results are clus- 
tered around M31's systemic recession velocity of —300 km 
S - . Regarding this issue, Chapman et al. (2006) analysed 
a large sample of red giant branch stars in the halo of M31 
(between 10 ^ R p ^ 70 kpc) and fitting a simple model to 
the observed data they obtained the dispersion of the radial 
velocity as a function of the projected radius -RpQ- Their lin- 
ear relation implies a VQ — 62 km s _1 at 100 kpc distance 
from the centre of M31. Combining the velocity results for 
the 5 globular clusters beyond R p = 70 kpc presented in 
Fig. [2] we derive a velocity dispersion of 65 km s _1 in the 
outer halo of M31. For the Milky Way globular cluster sys- 
tem, Battaglia et al. (2005) found that the radial velocity 
dispersion shows a nearly constant value of 120 km s" 1 out 

1 o V q(R p ) = 152 - 0.90i?,p km s , with R p given in kpc. 



to 30 kpc, which declines down to 50 km s" 1 at about 120 
kpc. The declining velocity dispersion in the outer halo sug- 
gests a relatively low mass halo, in which the virial velocity 
of the halo is less than the rotation speed of the disk (Abadi, 
Navarro & Steinmetz 2006), similar to that proposed for the 
Milky Way by Klypin, Zhao & Somerville (2002). 

In Fig.|3]we show the measured globular cluster internal 
velocity dispersion versus the absolute visual magnitude. It 
is evident that our 3 outer globular clusters resemble those 
found for Galactic globular clusters (Pryor & Maylan 1993) 
and also for other globular clusters in M31 (Barmby et al. 
2007). Such a scaling relation is not well understood to date. 
Yet, as briefly discussed by Djorgovski (1991, 1993), the pri- 
mordial L-a relation (slope ~1) of young globular clusters 
would be affected by subsequent dynamical processes (tidal 
shocks, for example), which would lead to the L-a power-law 
relation observed by the present day. 

For our sample, the dynamical mass-to-light ratio in 
the V band M/Ly ranges from 1.50 to 4.51 (in solar units). 
Our results are in agreement with those obtained by Du- 
bath & Grillmair (1997) for nine M31 globular clusters (1.0 
s: M/Lv ^ 6.5). In addition, Pryor & Meylan (1993) ob- 
tained a mean dynamical M/Lv = 2.3 ± 1.1 for the Galac- 
tic globular cluster system. Hence, within the uncertainties, 
GCM06 and GC5 display typical Galactic globular cluster 
M/Lv ratios, while GC10 is rather on the high side. 

When compared to the Galactic globular cluster sys- 
tem, one remarkable aspect of the M31 globular cluster sys- 
tem is its large range in age. That is, while the Milky Way 
presents a very homogeneous old population, M31 displays a 
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Table 3. Metallicity determinations for the 3 outer globular clus- 
ters. 



-10 



-12 



M v [mag] 



Figure 3. Internal velocity dispersion as a function of the abso- 
lute visual magnitude for the present sample (triangles) , Galactic 
globular clusters (open circles) and for other globular clusters in 
M31 (solid circles). Refer to the text for references. 



Table 2. Derived parameters for the 3 outer globular clusters. 



ID 


Rp 


(7 


Al 


M/L v 


Age 




[kpc] 


[kms — 1 ] 


[1O 5 M ] 


[Mq/Lq] 


[Gyr] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


GCM06 


118.2 


7.7 ± 0.4 


3.14 ± 1.57 


1.50 ± 0.75 


7.1 ± 3.0 


GC5 


79.4 


7.2 ± 0.4 


7.51 ± 1.50 


2.73 ± 0.54 


10.0 ± 3.0 


GC10 


101.7 


8.8 ± 0.7 


7.83 ± 1.56 


4.51 ± 0.90 


12.6 ± 3.0 



young population 2 Gyr), a number of intermediate age 
globular clusters (3-6 Gyr) and old globular clusters 7 
Gyr) as well (e.g. Beasley et al. 2005). Based on our stellar 
population analysis, we found that the present outer halo 
globular clusters of M31 have a mean age of around 10 ± 3 
Gyr. This mean spectroscopic age determination is in agree- 
ment with the color-magnitude diagrams presented in Mar- 
tin et al. (2006) and Mackey et al. (2007) for these globu- 
lar clusters. Furthermore, our dynamical mass-to-light ratios 
(M/Ly ^ 1.5) are compatible with an old stellar population. 
Our final derived parameters are summarized in Table [2] 

Spectroscopic metallicities for outer metal-poor globu- 
lar clusters are essential to the scenarios of chemical enrich- 
ment of the halo in galaxies since these objects are consid- 
ered halo-tracers. In Table [3] we present our spectroscopic 
metallicity determinations and, for comparison, the photo- 
metric metallicities. The final metallicities were computed 
as the median of the different spectroscopic values presented 
in Table |31 with a robust uncertainty of 1.5xSIQR (Semi- 
InterQuartile Range) being adopted. This lead to [Fe/H] = 
-1.37 ± 0.15 dexfor GCM06, [Fe/H] = -1.33 ± 0.12 dex for 



Method 


GCM06 


GC5 


GC10 


x 2 


-1.65 ± 0.10 


-1.40 ± 0.10 


-1.73 ± 0.10 


Mgb*> a 


-1.56 ± 0.11 


-1.74 ± 0.08 


-1.82 ± 0.08 


CH*' 6 


-1.37 ± 0.25 


-1.24 ± 0.25 


-1.81 ± 0.24 


Mg2*- C 


-1.36 ± 0.05 


-1.33 ± 0.06 


-1.54 ± 0.01 


CMD d 


-1.30 ± 0.15 


-1.84 ± 0.15 


-2.14 ± 0.15 


Adopted 


-1.37 ± 0.15 


-1.33 ± 0.12 


-1.73 ± 0.20 



Notes. — (*): Average values using different calibrations 
as follows : (a),(b): Brodie & Huchra (1990) and Perrett 
et al. (2002); (c): Brodie & Huchra (1990) and Buzzoni 
et al. (1992); (d): Isochrone fitting from Martin et al. 
(2006) and Mackey et al. (2007). 



GC5 and [Fe/H] = -1.73 ± 0.20 dex for GC10. These values, 
taking all uncertainties into account, are consistent with the 
CMD ones reported in Table|3]for GC10 and GCM06, while 
GC5 is a factor of 3 more metal-rich than the CMD esti- 
mate. Note, however, that although the metallicity quoted 
for GC5 based upon the Mgb method is also supported by 
visual inspection of the high-quality spectra shown in Fig. 
[T] the higher value obtained by the median of the different 
spectroscopic measurements was adopted. 

Both the photometric and spectroscopic methods have 
their advantages and disadvantages. Whilst the former is 
reddening and distance- modulus dependent, the latter is 
calibration and model-dependent. Furthermore, photomet- 
ric metallicities are also susceptible to the set of isochrones 
used as well as to age-metallicity degeneracy on the red giant 
branch. On the other hand, measurements of age, metallic- 
ity and elemental abundance from integrated spectra using 
a single index or pairs of orthogonal indices (e.g. H/3 vs. 
[MgFe]) are especially susceptible to index calibration er- 
rors, a particular concern for our uncalibrated data. 

In Fig. |3] we present the behaviour of the mean metal- 
licity as a function of the projected radii in the halo of M31 
using our metal-poor globular cluster sample, together with 
additional values from the literature for field stars and glob- 
ular clusters. Three main aspects are illustrated in this fig- 
ure. Firstly, there seems to exist a metallicity gradient for 
the metal-poor globular clusters at a projected distance less 
than 50 kpc from the centre of M31. Secondly, we find ten- 
tative evidence for the mean globular cluster metallicity to 
flatten off at [Fe/H] ~ -1.6, which is similar to the lower 
threshold seen in many other galaxies (e.g. Forbes et al. 
2000). Thirdly, our results confirm that the globular clus- 
ters in M31 are systematically more metal-poor than their 
outermost counterpart field stars, bearing in mind that we 
may be sampling different physical radii. This point can be 
explained in the context of chemical pre-enrichment hav- 
ing taken place in M31 before most of the halo formation 
has occurred (Durrell, Harris & Pritchet 1994, 2001). Simi- 
lar results have also been found for giant elliptical galaxies 
(e.g. Harris, Harris & Poole 1999). However, Schorck et al. 
(2008) show that the metallicity distribution function of the 
Galactic halo field stars is indistinguishable from that of the 
Galactic globular cluster system. 
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Currently, it is not clear what fraction of the Milky 
Way's halo has been assembled via dissipational collapse 
(monolithic collapse and gas rich mergers), and which by 
dissipationless accretion. However, the global properties of 
the present-day globular clusters of the Milky Way suggest 
a scenario where the inner halo (within 10 kpc) was formed 
via dissipative collapse, while the outer halo has an accretion 
origin (e.g. Bica et al. 2006). The metallicity trend in Fig. 
[4] is also consistent with this picture. It suggests that the 
steep metallicity gradient for the globular cluster system of 
M31 in the inner halo flattens in the outer regions. Such 
a flattening could be indicative of the outer halo globular 
clusters being likely accreted from satellites during an early 
period of merging (see e.g. Font et al. 2006). 



5 SUMMARY 

In summary, we confirm that the 3 outer globular clusters in 
the present sample are kinematically members of the M31 
globular cluster system and have structural parameters that, 
within the uncertainties, resemble those of the Milky Way 
globular clusters. Furthermore, they are characterized by an 
old stellar population (> 7 Gyr), whose mean metallicity 
matches the M31 metallicity peak for halo globular clusters 
(e.g. Barmby, Huchra & Brodie 2001). Unlike the Milky Way, 
we confirm that the outer halo globular clusters of M31 are 
systematically more metal-poor than their counterpart field 
stars (Kalirai et al. 2006) at a given projected radius. We 
also find evidence for a declining velocity dispersion with 
radius for the globular cluster system. 

We note that the most distant globular cluster in the 
Milky Way, AM 1, lies at 120 kpc. It is reasonable to expect 
that more globular clusters in the outer halo of M31 will 
be discovered with time, which should be followed-up with 
integrated-light spectra. 
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